EPJ manuscript No. 

(will be inserted by the editor) 



o : 
o . 

^ Magnetic Anisotropy of Deposited Transition Metal Clusters 

[i, ■ S. Bornemann 1 , J. Minar 1 , J. B. Staunton 2 , J. Honolka 3 , A.Enders 3 , K. Kern 3 and H. Ebert 1 



(N 



CZ3 



o 
o 



> 
oo 

(N 
O 

o 

-I— > 
c3 



i 

c 

o 
o 



X 



1 Department Chemie und Biochemie, Ludwig-Maximilians-Universitat Miinchen, 81377 Miinchen, Germany 

2 Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 

3 Max Planck Institut for Solid State Research, 70569 Stuttgart, Germany 

Received: / Revised version: 

Abstract. We present results of magnetic torque calculations using the fully relativistic spin-polarized 
Korringa-Kohn-Rostoker approach applied to small Co and Fe clusters deposited on the Pt(lll) surface. 
From the magnetic torque one can derive among others the magnetic anisotropy energy (MAE). It was 
found that this approach is numerically much more stable and also computationally less demanding than 
using the magnetic force theorem that allows to calculate the MAE directly. Although structural relaxation 
effects were not included our results correspond reasonably well to recent experimental data. 
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1 Introduction 

In recent years, low-dimensional magnetic nanostructures 
on surfaces have been the subject of intense experimen- 
tal and theoretical research activities which are driven 
by fundamental as well as practical interests. One of the 
central questions for potential future applications is how 
the spin-orbit coupling (SOC) induces specific magnetic 
properties such as orbital magnetic moments and mag- 
netic anisotropy. Another interesting issue is its influence 
on the exchange coupling. Hereby, it is especially impor- 
tant to understand how these properties evolve from sin- 
gle magnetic adatoms to submonolayer magnetic particles. 
This very prominent role of SOC for such systems is also 
reflected in recent reviews pQ as well as theoretical work 
0. 



theorem or by determining the total energies as a func- 
tion of the magnetisation direction. However, it turned 
out that these procedures are rather delicate when deal- 
ing with deposited clusters and one needs to take great 
care when taking band or total energy differences. There- 
fore, we implemented a method to calculate the magnetic 
torque directly from the electronic structure. Calculating 
the magnetic torque for a sequence of directions of the 
magnetisation then gives access to the MAE. 

In this present work we show first results for magnetic 
torque calculations of the already well investigated Co 
adatoms and dimers deposited on Pt(lll) and compare 
them with their Fe analogues. These investigations are 
complemented by calculations for decorated clusters that 
allow to optimize the MAE. Our theoretical results are 
then used to simulate magnetisation curves of an ensemble 
of Fc n (n=l,2,3) clusters on Pt(lll), that are compared 



A system which has been intensively investigated ex 
pcrimentally [3l4l5l6l7l8l9ll0] as well as theoretically [llllfllSFLlffSr * experimental results 
in recent years is Co/Pt(lll) (used here as a short nota- 
tion for Co clusters or nanostructures, respectively, de- 
posited on a Pt(lll) substrate) as this is a prototype 
' to study the requirements on new high-density magnetic 
storage materials. Earlier theoretical works studied only The calculations for the investigated cluster systems were 
rather small Co clusters or Co chains on Pt(lll) |llll2ll3U4] done within the framework of spin density functional the- 



2 Computational Details 



whereas only recently first qualitative results and trends 
based on a parameterised tight-binding approach were 
published for deposited structures of up to 37 Co atoms 
[T5] . In our earlier works [Illl2ll3j on the Co/Pt(lll) 
system we already described the evolution of the spin and 
orbital moments as well as the exchange coupling for small 
clusters. What has been missing so far, for a complete pic- 
ture of the magnetic behaviour is the magnetic anisotropy 
energy (MAE) for such systems. In principle, calculations 
of the MAE are possible by applying the magnetic force 



ory using the local spin density approximation (LSDA) 
with the parameterization given by Vosko, Wilk and Nu- 
sair for the exchange and correlation potential [IB]. The 
electronic structure is determined in a fully relativistic way 
on the basis of the Dirac equation for spin-polarised po- 
tentials which is solved using the Korringa-Kohn-Rostoker 
(KKR) multiple scattering formalism [17ll8j . This proce- 
dure consists of two steps. First the Pt(lll) host surface 
is calculated self-consistently with the tight binding (TB) 
version of the KKR method using layers of empty sites 
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to represent the vacuum region. This step is then followed 
by treating the deposited clusters as a perturbation to the 
clean surface with the Green's function for the new sys- 
tem being obtained by solving the corresponding Dyson 
equation. This scheme is described in more detail in earlier 
publications [12113] . 

In all calculations the cluster atoms were assumed to 
occupy ideal lattice sites in the first vacuum layer and 
no effects of structure relaxation were included. There- 
fore, our results contain a systematic error and are strictly 
spoken not directly comparable with experimental data. 
Nevertheless, it could be shown already in earlier work 
[11; 13 14 on deposited clusters that this approach is ca- 
pable of reproducing systematic trends as well as achieving 
a reasonable quantitative agreement with values found in 
experiment. 

The aj}-component of the torque vector T"i"' on the 
magnetic moments oriented along the direction n was cal- 
culated with an expression based on Lloyd's formula and 
perturbation theory [19] 

T Lf = ~ J dE fME) (1) 
Im^tr(ri ft) [(« • J&W -1 - • I)]) ■ 

i 

Here fpD is the Fermi distribution function, u is the 
direction of the torque vector and J is the total angular 
momentum operator. Finally, the matrices and r^" 
are the single site t-matrix and the site diagonal scatter- 
ing path operator, respectively, where i (used as an in- 
dex) labels the atomic sites. As the calculations are done 
assuming the temperature T=0K the Fermi distribution 
function is replaced by the theta function 0(Ep — E) with 
Ep being the Fermi energy. For all results shown below, 
the energy integral / dE occuring in Eq. (TTJ) was calcu- 
lated on a rectangular complex energy mesh containing 
64 points, while using an angular momentum expansion 
up to i max = 2 for all the occuring matrices. 

Eq. fl} uses the analytic derivative of the energy with 
respect to a rotation angle. We found that this approach 
is numerically much more robust than taking the differ- 
ences between band or total energies. The disadvantage, 
however, is that the magnetic anisotropy energy, defined 
as the difference E(n,ho) of the energy for two orienta- 
tions of the magnetisation, n and no, respectively, has to 
be determined by a corresponding path integral: 

E(n,n Q ) = / T [fl) dn . (2) 

J h 

Developing E(fi,ho) in spherical harmonics up to sec- 
ond order and taking into account the symmetry of the 
investigated cluster substrate systems, one finds for ex- 
ample for a cluster having C2 V -symmetry with respect to 
its spatial structure, i.e. ignoring the orientation n of the 
magnetisation [20] : 




Fig. 1. Structures of the investigated systems: The cluster 
atoms occupy ideal lattice sites of the underlying Pt(lll) sub- 
strate. 



E(9, cf)) = Eq + K% t \ cos 26 (3) 
+-^2,2(1 — cos 29) cos 2cf> + K2,3 sin 29 sin <j> . 

Using a corresponding expression for the torque it is 
straight forward to deal with the integral occurring in 
Eq. ([2]). The evaluation of the anisotropy constants K n ^ m 
occurring in this equation can then be determined in a 
rather easy way by determining the torque for certain 
orentations h i.e. at angles (9, </>) of the magnetic moments 
(see below). 

3 Results and Discussion 

The structure of the investigated Co and Fe monomers and 
dimers are shown together with the underlying Pt(lll) 
substrate in Fig. [TJ As the ad- or cluster atoms, respec- 
tively occupy regular lattice sites correspondingly to the 
substrate the resulting cluster/substrate system has C 3v - 
and C2v-symmetry, respectively. Comparing the resulting 
spin and orbital magnetic moments of Co and Fe in Tab.[T] 
one notices that the spin magnetic moments for Fe are in 
general about 1.5 times larger than for Co. For both tran- 
sition metals the dimer formation has only a minor im- 
pact on their spin magnetic moments when compared to 
the single adatoms. The orbital moments, however, show 
a more interesting behaviour. Here we find already a sub- 
statial quenching when going from single adatoms to the 
corresponding pure i.e. unmixed dimers. This effect seems 
to be much more pronounced in the case of Fe where the 
orbital magnetic moment of a Fe2 dimer atom reduces to 
about one third of the monomer value compared to only 
three quarters in the case of Co. The corresponding val- 
ues for the mixed dimer differ only slightly from those of 
C02 and Fe2, respectively. The different behaviour of the 
orbital moment for Co and Fe is also reflected in their 
anisotropics (see below). A further increase in cluster size 
leads usually (depending also on the cluster shape) to a 
rapid and monotonous decay of spin and orbital magnetic 
moments approaching rapidly the values of a monolayer. 

Fig. [2] shows the dependence of the ^-component of 
the magnetic torque Tg(9, 4>) = — ^ on the polar angle <j> 
for Co and Fe monomers at 9 — j. As Tg(9,4>) is found 
to be negative here as well as for the following we show 
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Table 1. Spin and orbital moments per atom of the investi- 
gated systems with magnetisation along the z-axis, i.e. perpen- 
dicular to the surface. For the mixed cluster the value refers to 
the component underlined. 





Cor 


Fei 


Co 2 


Fe 2 


FeCo 


FeCo 


A'spin (Mb) 
Morb (/*b) 


2.27 
0.60 


3.49 
0.77 


2.15 
0.44 


3.33 
0.24 


2.12 
0.39 


3.38 
0.26 




Fig. 2. Torque component Tg(9, <j>) for Coi (top) and Fei (bot- 
tom) at 9 — j as a function of the polar angle <j> (In fact 
—Tg(6,<j>) = is shown - see text). The dashed line shows 
the same results including the substrate effect. 



always ^ = —Tg(9,(f>). The positive sign of ^ for all an- 
gles <p implies that the torque forces the magnetisation to 
the z-axis. This means that the system's easy axis points 
out of plane along z. One can see that the threefold sym- 
metry imposed by the underlying Pt substrate is directly 
reflected by the small oscillations of Tg with the polar an- 
gle <f>. 

This also demonstrates the high sensitivity of our im- 
plementation. The curves show no numerical noise even 
for an energy resolution below 0.1 meV. As the <j> depen- 
dence of Tg is so small when compared to its absolute value 
the adatoms almost behave like perfect uniaxial magnets. 
In this case the anisotropy constants K can be extracted 
from the minima of T$(9,(f)). For the single adatoms this 
gives then 4.76 meV for Co and 8.79 meV for Fe. Fig. [5] 
also shows the influence of the induced anisotropy coming 
from the Pt substrate atoms. This induced MAE is about 
-30 /xeV for Fe and even smaller in the case of Co. In fact it 
seems to be negligible for very small Co clusters composed 




Fig. 3. Top panel: torque component To (9, 4>) for C02 (circles), 
Fe2 (squares) and the FeCo dimer (thick line) for 9 = ? . Lower 
panel: (f> dependence of the average orbital moments for all 
three dimers at 9 = ^ . 



of only few atoms. For larger two-dimensional Co clusters, 
however, we found that this induced anisotropy becomes 
more important with increasing cluster size and can rise to 
the same order of magnitude as the contribution coming 
from the Co atoms themselves. 

The results of the torque calculations for the three dif- 
ferent investigated dimers are shown in the top panel of 
Fig- El In all cases ^ is again positive, however, the abso- 
lute values are significantly reduced when compared to the 
monomers. The (^-dependence of the ^-component of the 
torque on the other hand is now much more pronounced 
due to the reduced symmetry of the cluster/substrate sys- 
tem when compared to the monomer (see Fig. [T]). From 
Fig. [3] one can see that E(h, i) is smallest if the mag- 
netic moments are oriented along the x direction, i. e. 
along the cluster dimer axis (see Fig.QJ. The lower panel 
of Fig. [3] shows the corresponding <fi dependence of the 
orbital moments for 6 = \. One can see that the oscilla- 
tions in the orbital moments follow the oscillations of Tg in 
an anticyclic manner. Here it should be pointed out that 
the largest values for the orbital moments are obtained 
when the magnetisation points along the z-axis (see Ta- 
ble [1]) . This as well as the behaviour seen in Fig. [3] is in 
qualitative agreement with the Bruno and van der Laan 
anisotropy models [21122) that relate the MAE E(h, no) 
to the corresponding anisotropy of the orbital moment. 

For the van der Laan model one has 
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Fig. 4. Top panel: Torque component T,f,(9, <j>) of C02 for 9 = 5 
as a fuction of 0. Middle panel: MAE AE(0) for 6 1 = f refered 
to no— (■§, = 0). The thick line gives calculated results while 
the thin line marked by diamonds gives results based on van der 
Laan's model. Lower Panel: spin resolved anisotropy Afi m3 ((f>) 
for 9 — ^ refered to no=(f,^> = 0) for spin up and down 
(triangles up and down, resp.) and summed up (diamonds). 



E(n,n ) = --^[^(n) - fxl h (n)] (4) 

4^B 

+ 7— [Morb(«o)-Mi rb ("o)] , 
4/i_B 

where £ is the appropriate spin-orbit coupling constant 
(here for the 3d transition metal Fe or Co, respectively), 
and ^™^(n) is the spin projected orbital magnetic moment 
for an orientation of the magnetisation along n. If one con- 
siders a strong ferromagnet one can ignore the majority 
spin contribution in Eq. (j4|) leading to the expectation 



Fig. 5. Structure of the studied Fe trimer on Pt(lll) without 
(left) and with (right) decoration by three Pt atoms. 

that the orbital magnetic moment takes its maximum for 
the magnetisation oriented along the easy axis. The results 
shown in Figs. [2] and [3] are obviously in full agreement with 
this. As Fig. [3] shows for the three considered dimers, also 
the MAE deduced from the shown torque goes parallel 
with the orbital magnetic moment (lower panel) when <f> 
is varied while 9 is kept fixed. 

The relation of the MAE and the anisotropy of the 
orbital magnetic moment has been studied in more detail 
by calculating the 0-component of the torque for 9 fixed 
to ^, i.e. for the magnetisation forced to lie in the surface 
xy-plane. As one notes from the results given in Fig. [4j 
the variation of the torque is again quite pronounced and 
reflects the C2 V -symmetry of the dimer/substrate system. 
Integrating the torque component with respect to 4> ( see 
Eq. ©) one is led to the MAE AE((f,c?!>), (§,</>)) shown 
in the middle panel of Fig. |H that again shows that the 
energy is at its minimum if the magnetisation lies along 
the dimer axis. The corresponding anisotropy Afi(<fi) is 
given in the lower panel of Fig. [4] in a spin-polarized way. 
As one notes the majority contribution to Z\/z or b(0) = 
A*orb(f ,<t>) — Morb(f ,0) is very small and negligible. This 
is a consequence of the nearly complete filling of the ma- 
jority spin d-band of the Co-atoms, i.e. their strong fer- 
romagnetic behaviour. As a consequence, one has not to 
distinguish between Bruno's and van der Laan's models. 
Using the relation given in Eq. ^ together with the spin 
orbit coupling strength £ = 85 meV for Co [23] one is 
lead to the estimate for the MAE AE(cf>) in the middle 
panel represented by the curve marked by diamonds. Ob- 
viously, the qualitative behaviour of the MAE is properly 
reproduced with a reasonable quantitative agreement. 

Figs. [2] and [3] show that going from a monomer to a 
dimer the MAE is drastically reduced. This goes on with 
increasing cluster size, although there is some influence of 
the cluster shape [3]. As one notes, the decrease of the 
MAE is much more pronounced for Fe than for Co. For 
the mixed FeCo cluster, however, the MAE is quite high. 
This indicates that by a suitable combination of atoms one 
may optimize the MAE while keeping the magnetisation 
high. In fact, there are already a number of experimental 
studies done along this line. Here, we show results for the 
impact on the properties of a Fe trimer on Pt(lll) due to 
a decoration with Pt atoms. 

Fig. [5] shows the corresponding structure of the inves- 
tigated cluster/substrate systems. The resulting spin and 
orbital magnetic moments are given in Table O As one 
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Table 2. Spin and orbital moments per atom of the Fe3 and 
FeaPt3 clusters on Pt(lll) for the magnetisation along the z- 
axis, i.e. perpendicular to the surface. The values refer to the 
component underlined. The last two columns give results for 
the Pt-substrate atoms having one or two Fe atoms as nearest 
neighbours, resp., in case of Fe3Pt3. 





Fe 3 


Fe 3 Pt 3 


Fe 3 Pt 3 


Pt(l) 


Pt(2) 


Atspin (Mb) 


3.20 


3.16 


0.21 


0.09 


0.11 


Morb (fJ-B) 


0.16 


0.09 


0.08 


0.02 


0.02 



notes, the moments of the Fe atoms are again quite high 
as for the mono- and dimer. For the spin moment one finds 
only a minor influence due to the decorating Pt atoms. 
These have also quite an appreciable induced magnetic 
moment. The resulting torque Tg(9, (f) for 9 = ? as a func- 
tion of the angle 4> is shown in Fig. [6] (top panel) for the 
Fe2-, Fe3- and Fe3Pt3 clusters. As is demonstrated once 
more the torque Tq and with this also the MAE is strongly 
reduced going from the dimer to the trimer. Adding the 
decoration, however, the high MAE of the dimer is re- 
covered. As the inner Fe3 cluster of the Fe3Pt3 cluster is 
now surrounded by neighbouring atoms the variation of 
the MAE with <fi is strongly reduced when compared with 
the dimer. 



As one would expect, the MAE of the decorated trimer 
Fe3Pt3 is dominated by its Fe contribution. This can be 
seen in the middle panel of Fig. [5] where apart from the 
total MAE per Fe-atom the averaged contribution of a Fe- 
atom is shown. Nevertheless, there is an appreciable con- 
tribution coming from the decorating Pt atoms as well. 
This has to be ascribed on the one-hand side to their ap- 
preciable induced magnetic moment (see Tab. [2]) and on 
the other side to their large spin-orbit coupling (£ Pt =710 
meV to be compared with £pc=65 meV). Apart from the 
contribution of the decorating Pt atoms there is a note- 
worthy contribution from the neighbouring Pt-substrate 
atoms as well. As can be seen in the lower panel of Fig. [6] 
this amounts to about 5% of the total torque or MAE, 
respectively. 

The results for the undecorated Fe„ clusters (n=l,2,3) 
shown in the lower panel of Fig. [2] and the top panel of 
Fig. [6] allows us to make contact to corresponding exper- 
imental investigations |24j . The sample preparation tech- 
nique used within these investigations led to an ensemble 
of Fe„ clusters that was dominated by clusters of size n=l- 
3 with their statistical weight determined to be w„=0.53, 
0.54 and 0.33 % (for technical details see Ref. 0). Fig. [7] 
shows the magnetisation curves M(B) measured for this 
cluster ensemble at T=6K for an orientation of the exter- 
nal magnetic field B along the easy axis (z) and at an angle 
9=65° with respect to this axis. With the magnetic mo- 
ments and the anisotropy parameters for the Fe n clusters 
available the magnetisation curves M(B) can be simulated 
by means of the so-called Langevin formula 4 . This way 
the thermal average of the z-componcnt m nz (B,T) of the 
moment m n of an Fe„ cluster can be expressed by: 




Fig. 6. Top panel: torque component Tg(8, cf>) for Fe2 (circles), 
Fe3 (triangles) and Pt3Fe3 (squares) for 8 — ? as a function 
of the polar angle (f>. Middle Panel: torque component Tg(8, <j>) 
per Fe atom for Pt3Fe3 at 8 = j as a function of <j>'- total (full 
squares) including the substrate contribution, contribution of 
a Fe atom (triangle) and a cluster Pt atom (circles). Lower 
Panel: total torque component Tg(8, (j>) per Fe atom for Pt3Fe3 
at 8 = j as a function of the polar angle </> with (full squares) 
and without (open squares) the substrate contribution. 



/* sin 9d9e- E ^ B ' T ' e y kT m n cos 

mnz{B ' T) ~ jf^m^wmikf • (5) 

For the simulation the energy E(B, T, 9) was assumed 
to consist of its Zeeman and anisotropy contributions 



E(B, T, 9) = Bm n cos 9 + sin 2 9 



(6) 



where for the later one an uniaxial behaviour has been 
assumed. The corresponding anisotropy constant K" has 
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Fig. 7. Experimental magnetisation curves M(B) (dots) of 
an ensemble if Fe n clusters on Pt(lll) measured at T=6K 
for an orientation of the magnetic field M(B) along the easy 
axis z(6 = 0°) and rotated by 6 = 65° with respect to that. 
The full lines give corresponding theoretical results obtained 
on the basis of the calculated properties of Fe„ clusters and the 
Langevin formula given in Eq. ([5]). The dashed line is obtained 
by including Fe4 clusters in the simulation. 

been deduced from the results for the MAE shown above. 
Adding the magnetisation curves for the Fe„ clusters weigh- 
ted by their statistical weight w„ one is lead to the full 
lines shown in Fig. [7] The additional dashed lines stem 
from a second simulation done including tetramers. This 
indicates that a certain amount of larger Fe clusters are 
formed during the preparation process as expected by 
statistics. 

The good agreement of the simulated curves with ex- 
periment obviously demonstrates that the complex mag- 
netic properties of transition metal clusters can indeed be 
understood and described without adjustable parameters 
on the basis of the approach used here. It also shows that 
inclusion of relaxation effects when calculating the cluster 
properties should only slightly modify the numerical re- 
sults. Nevertheless, corresponding numerical works are in 
progress to determine the influence of lattice relaxations. 

4 Conclusion 

We have presented results for the magnetic anisotropy of 
various Co/Pt(lll) and Fe/Pt(lll) cluster/substrate sys- 
tems calculated by the fully relativistic KKR-approach. It 
was demonstrated that calculating the magnetic torque 
instead of the magnetic anisotropy energy directly, is nu- 
merical very robust and reliable. This is in particular re- 
flected by the accuracy achieved for the dependency of the 
torque on the polar angle 4> as well as the results for the 
substrate contribution. 

As in previous work it was found that increasing the 
cluster size in general a rapid decrease of the MAE occurs. 
However, it could be shown that by suitable formation 
of compound clusters this drop can be compensated. In 
fact, the combination of an element with large magnetic 



moments with one having large spin-orbit coupling seems 
to be a promising approach. 

Finally, it could be demonstrated that using the calcu- 
lated cluster properties the results of experimental mag- 
netisation curves could be reproduced in a very satisfying 
way confirming the adequateness of our approach as well 
as the interpretation of the experimental findings. 
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